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We analyze the moments of parton distribution functions in the pion calculated in lattice QCD, 
paying particular attention to their chiral extrapolation. Using the lowest three non-trivial mo- 
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ments calculated on the lattice, we assess the accuracy with which the x dependence of both the 



valence and sea quark distributions in the pion can be extracted. The resulting valence quark 
distributions at the physical pion mass are in fair agreement with existing Drell-Yan data, but the 
■ statistical errors are such that one cannot yet confirm (or rule out) the large- x behavior expected 



from hadron helicity conservation in perturbative QCD. However, one can expect that the next 
generation of calculations in lattice QCD will allow 
of accuracy comparable with current experiments. 



generation of calculations in lattice QCD will allow one to extract parton distributions with a level 
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I. INTRODUCTION 



It is widely appreciated that the pion plays a very fundamental role in QCD. Given that 
chiral symmetry is such a good symmetry of nature, because of the extremely low masses 
of the u and d quarks, the pseudo-Goldstone character of the pion is ubiquitous in hadron 
physics. As a result, the determination of its structure, both from experiment and non- 
perturbative studies of QCD, is of great importance. The parton distribution functions 
(PDFs) of the pion have been measured in a number of experiments, using the Drell-Yan 
reaction 0, 0, U, 0]. Such experiments tend to focus on the region of Bjorken-x above 
pa 0.2 and hence are most sensitive to the valence distribution. Until recently there was 
little constraint on the size or form of the sea quark distributions, but measurements of 
charge-exchange in semi-inclusive deejo inelastic scattering (DIS) at HERA have yielded 
some information at very low x |6J, 17J, |8| , and one can also expect new, high precision data 
from semi-inclusive DIS after the upgrade at Jefferson Lab 9]. This observation will be 
important for our analysis because the current errors for the sea quark distributions are 
considerably larger than the statistical errors in the first moment of the lattice data. 

The existing data have been used to constrain various phenomenological parameteriza- 
tions of the pion PDFs . At the same time they are used to guide and 

test non-perturbative models of the internal structure of the pion, from the constituent 
quark model Q to the NJL model Q, H H Q and others 
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there has recently been a calculation within a covariant model, based on a truncation of the 
Dyson-Schwinger equations 

One of the clearest predictions for the x dependence of the pion structure 
from considerations of hadron helicity conservation within perturbative QCD 
It is a firm expectation within this framework that the valence quark distribution should 
behave like (1 — x) 2 as x — > 1. On the other hand, the experimental data seem to be more 
consistent with a form linear in (1 — x). One suggestion is that the experimental data may 
have a substantial higher-twist component j^. We shall see that the analysis of data from 
lattice QCD offers a significant possibility of resolving the issue in the near future. 

In Section |H] we review the lattice simulations of the moments of the pion structure 
function, while the chiral extrapolation of these moments is described in Section IIHI The 
reconstruction of the x dependence of the valence and sea quark distribution functions in 
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the pion is presented in Section IIVI In order to make quite clear what can be learned from 
existing lattice data, and what might become possible in the near future, we present several 
alternative methods for performing the extraction. In Section IV, we also investigate the 
pion mass dependence of the reconstructed distribution. Finally, in Section |V] we summarize 
our results, and outline future applications of the methodology presented here. 



II. LATTICE RESULTS 



By discretizing space-time as a four- dimensional hyper-cubic lattice, the field equations 
of QCD can be solved numerically in the non-perturbative region. The potential of lattice 
QCD is that it allows a first principles investigation of hadron properties and structure. The 
main weakness of these numerical calculations is the vast computational resources that they 
require. Indeed, it is not yet computationally feasible to perform lattice calculations that 
correspond to the parameters of the real world. Current simulations are run at quark masses 
3-10 times too large, on lattice volumes that are likely too small, and often use the quenched 
approximation (in which sea quark loops are neglected). The result of these restrictions is 
that various extrapolations are necessary to reach the physical regime. 

The pioneering lattice calculations of hadron structure functions were made by Martinelli 
and Sachrajda in the late 1980s js^l^l- Even though the available computational resources 
restricted the statistical accuracy of their studies and confined them to small lattices, their 
results are still consistent with the more advanced calculations of the QCDSF collaboration 
which we discuss below. First, however, we briefly consider the formalism needed to connect 
the lattice and continuum theories. 

While the x dependence of the parton distribution functions cannot be computed directly 
on the lattice, one can compute the moments, (x n ) , of the distributions. Using the operator 
product expansion, these moments can be related to matrix elements of operators of a given 
twist. The leading twist (twist-2) operators are given by 

W-Mn = jn-1 ^ 7 {Mi D m . . . D^lffg , (1) 

where ip q are quark fields, _D M is the covariant derivative, and the braces {■ ■ •} denote sym- 
metrization of indices. For reference, we shall work with the u quark distribution in the ir + 
meson, (x) , which can be related to distributions in the 7r~ and 7r° by charge symmetry 
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(c.f. Refs. (32,|33(), 



u w +(x) = d w +(x) 



d w -(x) 



q n {x) 



(2) 



where we have suppressed the dependence on the scale Q 2 . The moments of q-n{x) are defined 
as 



where, for example, the n = moment corresponds to the number sum rule, (x°) q = 1. 
Operationally, these moments can be extracted from the forward matrix elements of the 
operators (JH) as 



where "traces" are subtracted to give matrix elements that transform irreducibly. 

In the lattice formulation, discretized versions of the operators (P) must be defined that 
have the correct continuum limits. A number of technical considerations arise in this pro- 
cedure. For the n = 1 moment, there are two possible lattice discretizations of the corre- 
sponding continuum operator. One of these can be evaluated with both pion states having 
zero momentum, which results in greater statistical precision. We only include the results 
for this operator in our analysis. The lattice data for the less well determined operator are 
consistent with this, however, and their inclusion would not modify our conclusions. For 
n = 2 and 3, non-zero momentum is unavoidable and the data are correspondingly less pre- 
cise. Also, the reduced symmetry of the lattice means that it becomes impossible to define 
operators that transform irreducibly for n > 4. Calculation of the corresponding moments 
is more difficult as it necessarily involves the evaluation of coefficients which describe the 
mixing with lower dimensional operators. Consequently there are only data for n = 1, 2 and 
3 at the present time. 

Although somewhat easier to calculate than for the case of the nucleon, the moments of 
the pion distribution functions have received less attention in the literature. The QCDSF 
collaboration has performed the only detailed study j^J of the moments of the pion parton 
distributions. The analysis was based on a sample of 500 configurations, with the simulations 
performed in the quenched approximation using a Wilson quark action at three different 
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quark masses, m q ~ 70, 130 and 190 MeV, on a 16 3 x 32 lattice at f3 = 6.0. QCDSF set 
the scale by linearly extrapolating the p meson mass to the chiral limit. Although there are 
considerable uncertainties associated with such an extrapolation, given the potential non- 
linearities associated with chiral non-analytic behavior, the study by Leinweber et al. 35^ 
suggests that a linear approximation may not be so inaccurate in this particular observable. 
In any case, with the physical scale set in this way the lattice moments correspond roughly 
to a scale Q 2 « (2.4 GeV) 2 ~ 5-6 GeV 2 Q- The QCDSF collaboration have also analyzed 
some higher twist contributions to the pion structure function 36| . finding that they are 
rather small (at least at the large quark masses considered). 

While the QCDSF investigations used quenched field configurations, one would expect 
that the effects of that approximation should be relatively small at the large quark masses 
for which data are available. Indeed, prev ious comparisons of quenched and unquenched 
data for nucleon structure calculations I37L showed no statistically significant difference 
in this region. The QCDSF lattice study of hadron structure is ongoing and we look forward 
to unquenched results in the near future. When lattice calculations are able to be performed 
at significantly lighter masses, the effects of quenching will become apparent. Finally, we 
note that the lattice results cannot be regarded as definitive, even at the masses used, until a 
thorough investigation of the effects of the finite lattice spacing and finite lattice volume has 
been undertaken. For example, Jansen [3^ suggests that the 0(a) errors could be significant 
in calculations of (x) q with Wilson fermions. Bearing these caveats in mind, we take the 
lattice results at face value in the current study, with the understanding that our analysis 
can easily be updated to reflect improvements in the lattice data as they occur. 

We stress that, even though lattice QCD calculations in the next few years will be ex- 
tended to smaller quark masses and larger lattices in the quenched and unquenched (or at 
least partially quenched 4(|) versions of QCD, the numerical challenge of light quark masses 
is such that extrapolation over a fairly large range of quark mass will be needed for many 
years. 
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III. CHIRAL EXTRAPOLATION 



The approximate chiral symmetry of QCD leads to the appearance of pseudoscalar (Gold- 
stone) bosons. In the case of chiral SU(2) i xSU(2) fl , these are identified with the pions, 7r ± '°. 
Because the pion mass vanishes with the square root of the current quark mass, m n ~ y/m q , 
the pion takes on an increasingly important role in QCD as m q — > 0. Its effect on hadron 
structure can be quantified using systematic expansions of observables in powers (and log- 
arithms) of m n |4l| . In particular, because of the structure of the Goldstone boson loop 
corrections to hadronic properties, coefficients of terms in the expansions which are non- 
analytic in the quark mass can be calculated in terms of physical parameters, and hence are 
model independent. For the case of the nucleon, the leading non-analytic behavior of the 
moments of parton distribution functions arising from such loops was found to be crucial 
in understanding the relation between the lattice results and the physical values of the mo- 
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47, 48]. lattice calculations from 

the unphysically large quark masses at which they are currently performed to the physical 
quark masses must incorporate the effects of the pion cloud [49]. 

Arndt and Savage j^J have calculated the leading chiral corrections to the moments of the 
pion's quark distributions, finding that the pion cloud contributions to the C-odd (n-even) 
flavor non-singlet (NS) moments receive corrections: 



n\NS 



(X»> 



q 



l-5 n0 2 ml 
;m„ log — 



(5) 



where f n = 93 MeV is the pion decay constant, a n is the value of the moment in the 
chiral limit and A x ~ ~ 1 GeV is the chiral scale. The n = moment is not 

renormalized by pion loops because of charge conservation. In the singlet sector, for the C- 
even (n-odd) moments, pion loops do not introduce any non-analytic structure. Physically, 
this is because any momentum lost by valence quarks through pion emission is recovered 
through the additional sea quarks generated. Of course, the C-even non-singlet and C-odd 
singlet moments must vanish identically because of the crossing symmetry properties of the 
distributions. 

Since the lattice data for the moments of the pion PDF are well fit by a linear function 
of m\, over the region where they have been calculated, it is natural to apply a functional 
form similar to that used to extrapolate the moments of the NS PDF of the nucleon j^. 
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We modify the linear term only minimally, replacing m\ by m^/(m^ + M 2 ) so that this term 
goes to a constant, rather than diverging, as — > oo, 



NS 



1 - c LN Aml log 



2 

mi 



y ml + /i 2 / 

2 



+ M z 



— 777 

where a n ,b n ,a n and 6 n are fit parameters, and clna = ^-/(^fn) 2 is the model independent 
coefficient of the leading non-analytic (LNA) term in the non-singlet expansion. The fits 
are insensitive to the parameter M as long as it is large, and in this analysis it is fixed at 
M = 5 GeV. 

The behavior of the moments in the limit m q — > oo can be determined model indepen- 
dently from heavy quark effective theory, so a more ambitious scheme would be to build this 
behavior into the fitting function as well. In the heavy quark regime, contributions from the 
quark- ant iquark sea are suppressed as and the two valence quarks in the pion each 

carry half of the momentum of the pion. The corresponding valence distribution is therefore 
a 5-function located at x = 1/2, so that the moments behave as 

(% n ) q — > 7^ , rn q -> oo . (8) 
This limit is easily built into the (non-linear) extrapolations, along with the chiral non- 



analytic behavior (see Ref. |43( for the analogous case of the nucleon). However, given the 
present accuracy of the lattice data it is sufficient to use the simpler extrapolation functions, 
given in Eqs. (jJJ) and (JBjl . which are not constrained by the heavy quark limit. 

The parameter \i in the argument of the chiral logarithm in Eq. (JBJ) is physically related 
to the size of the source of the pion cloud and controls the onset of the chiral behavior in 
the NS moments as m n —>■ 0. Ideally its value will be determined from fits to unquenched 
lattice data, however, present data are not yet at sufficiently low masses. Instead we take 
the value fx = 0.7 GeV, which is somewhat larger than that used in the nucleon analysis 
because of the smaller size of the pion, and test the sensitivity to fi by varying it over the 
range (0.4, 1.0) GeV. 

The above results for the chiral extrapolation are valid in full QCD, whereas the existing 
lattice data have been generated within the quenched approximation (in which the effects 
of background quark loops are neglected). Because quark loop effects are proportional to 
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FIG. 1: Chiral extrapolation of the lowest three lattice moments |2J] of the pion distributions. 
The upper plot shows the extrapolation of the valence moments and lower plot shows that of the 
total moments. The solid curves correspond to fits using Eq. ©, with ^ = 0.7 GeV for the valence 
moments, and Eq. for the singlet moments. The dark shaded region in both plots corresponds 
to fits to the data plus or minus their error bars, while the lighter shaded regions in the valence plot 
show the additional effect of varying ^ between 0.4 and 1.0 GeV on top of the statistical variation. 
The phenomenological valence (open stars) and total (open triangles) moments are shown at mP hys 
(see Section HVj) . g 



(powers of) l/m q , one expects loops to play a relatively minor role at large quark mass. 
Indeed, for moments of the nucleon parton distributions the quenched and full QCD sim- 



ulations were found 38( to be equivalent within statistical errors for m n > 0.5-0.6 GeV. 
Therefore, in the present analysis we assume that the available (quenched) data at large 
provide a reliable estimate of the unquenched moments at > 0.5-0.6 GeV. Future 
simulations will allow quantitative tests of this assumption, and when they can be performed 
at quark masses light enough for the difference to become apparent our analysis will need 



to be repeated using quenched [50j] (or partially quenched [5l() chiral perturbation theory. 

In general the matrix elements receive contributions from diagrams in which the operator 
insertions are either on quark lines which are connected to the pion source (CI), or on 
quark lines which are disconnected (DI) (i.e. connected only through gluon lines to the pion 
source) : 

(x n ) q = (xX + (x n )f ■ (9) 

The disconnected insertions contribute only to the singlet operators, while the connected 
diagrams contribute in both the singlet and non-singlet cases. The evaluation of disconnected 
diagrams is considerably more difficult numerically, and thus far only the connected pieces, 
(x n )^ 1 , have been computed [34J]. Once again we can make use of the large quark masses at 
which the lattice moments have been simulated by noting that the disconnected insertions 
should also be suppressed for ro, > 0.5 GeV, so that 

(x n ) q w (x n )f , m n > 0.5 GeV . (10) 

This same argument also suggests that, at these values ofm^, the pion PDFs should satisfy 

qtot a \(x) = q n (x) + q n (x) w q n (x) - q n (x) = v n (x) , (11) 

where v n (x) is the valence quark distribution in the pion. This observation allows one 
to approximate the C-odd (valence) moments by the C-even moments at large m^, and 
extrapolate them according to Eq. (jUJ) to compare with the phenomenological valence (non- 
singlet) moments. However, we stress that future lattice data for the n-odd moments should 
vary smoothly as m q decreases - i.e. they should show no rapid, non-analytic behavior as 
the chiral limit is approached. 

In Fig. [Jwe show the lattice data from the QCDSF collaboration for the n = 1, 2 
and 3 moments (of the u quark distribution in the tt + ) as a function of m^. The fits to 



the n = 2 data using Eq. © and those to the n = 1 and 3 moments using Eq. (J7J) are 
indicated by the curves in the upper and lower plots, respectively. For each fit the dark 
shaded error bands correspond to fits to the data ± errors. The phenomenological values 
of the moments, indicated by open stars (valence distribution) and open triangles (total 
distribution) at m^ hjs , are taken from an average of global fits 
function data (see Section Ft V Al below) . Assuming valence quark dominance of the moments 
at m n ^ 0.5 GeV, we also show in the upper plot the n = 1 and 3 moments extrapolated 
as if they were non-singlet, using Eq. ©• Under the same assumption we extrapolate the 
n = 2 moment as if it were a singlet in the lower plot. In the central curves of the NS fits, we 
choose fi = 0.7 GeV. The outer, lightly-shaded envelopes show a conservative variation of 
this parameter between 0.4 and 1.0 GeV in addition to the statistical variation (dark shaded 
region). In all cases the extrapolated moments, both singlet and non-singlet, agree with 
the phenomenological values within errors, as shown in Table H] This provides a posteriori 
evidence for the valence dominance of the moments (suppression of quark loops) at large 




to the pion structure 



IV. RECONSTRUCTION OF THE QUARK DISTRIBUTION 



In this section we use the available lattice moment data to constrain the Bjorken-x de- 
pendence of the underlying PDFs. The approach adopted here is similar to that in the 



earlier analysis of the PDFs in the nucleon 



44, 



521 ] . The general procedure is to choose a 



particular parameterization for the x dependence of the distribution, and perform a Mellin 
transformation to give the parametric dependence of its moments. Values for the moments, 
extrapolated from the lattice data, can then used to fit the various parameters and recon- 
struct the physical distribution. 



A. Phenomenological distributions 

Before using a specific parameterization to analyze the lattice data, we first test the 
robustness of the procedure by examining the extent to which the parameters of the phe- 
nomenological valence distributions can be reconstructed from their lowest moments. This 
will provide an estimate of the systematic error in the choice of parameterization and the 
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n = n = 1 


n = 2 


n = 3 


Moments of Phenomenological PDFs 


valence 


1 0.21(2) 


0.09(1) 


0.052(5) 


sea [Eq. ifTBJl] 


0.05(3) 


0.007(4) 


0.002(1) 


total 


0.31(6) 


0.11(1) 


0.056(6) 


Extrapolated Lattice Moments 


valence [method (ii)] 


1 0.24(1)(2) 


0.09(3)(1) 


0.043(15)(3) 


valence [method (iii)] 


1 0.18(6) 


0.10(3)(1) 


0.05(2) 


sea 


0.03(1) 




0.001(9) 


total 


0.275(8) 


0.11(3) 


0.05(2) 



TABLE I: Moments of PDFs of the pion, obtained from phenomenological PDFs (see Section HVj) 
and extrapolated from the lattice (as discussed in the text). The n = 2 lattice total moment is 
obtained from the lattice valence moment by adding twice the phenomenological sea. The lattice 
sea is determined by subtracting the phenomenological valence moments from the lattice total 
moments (for n-odd). Errors on the extrapolated lattice moments are calculated from fits to the 
data ± their errors (first parentheses) and from varying the parameter /i between 0.4 and 1.0 GeV 
(second parentheses, where applicable). 

reconstruction procedure. 

Several groups have performed global analyses of pion structure function data and con- 



structed parameterizations o 
the SMRS parameterization 



the parton distributions. The valence quark distribution in 
llT ] is fitted with the form 

v*(x) = Ax\l -x) c , (12) 



with the parameters A, b and c determined at an input scale of Q 2 = 4 GeV 2 , given in 
Table |H] According to Regge theory, the parameter b, which controls the x — > behavior, 
is given by the intercept of the p meson trajectory, and is predicted to be b ~ —1/2. The 
parameter c dictates the asymptotic behavior as x — > 1, and is pre dicted by hadron helicity 
conservation in perturbative QCD to be c = 2 for the pion The GRS (next-to- 



ll 



leading order) parameterization [12| contains two additional parameters, 

v n (x) = Ax b (l -x) c (l + e v / x + gx) , (13) 

with all parameters listed for Q 2 = 5 GeV 2 in Table |HI The small differences in scale 
between the parameterizations and the lattice moments are negligible. 

The phenomenological valence moments with which the lattice calculations are compared 
are defined by averaging the integrals of these two distributions, and the errors are calculated 
as the difference between the moments of the two distributions. These average moments are 
given in Table H] and shown at the physical pion mass as open stars (valence) and open 
triangles (total) in Fig. [T] 

The Mellin transform of the (more general) parameterization in Eq. (J13|) is given by 

(x n ) val = A[(3(l + c,l + b + n) + e (3(1 + c,3/2 + b + n) 

+g (3(l + c,2 + b + n)] , (14) 

where (3 (a, b) is the /^-function. For the simpler SMRS parameterization, only the first term 
in Eq. ()14j) is present. 

To determine our ability to reconstruct the parameters of a distribution from its moments, 
we first calculate the moments of the GRS distribution (to be specific) by direct numerical 
integration. Using Eq. (J14)) . we find that the five parameters in Eq. (|14j) can be very 
accurately reconstructed (to 4 significant figures) from the first five moments (n = — 4) 
using a standard Levenberg-Marquardt non-linear fit. However, since only 3 non-trivial 
lattice moments are currently available, one cannot determine all of the five parameters 
from the lattice data. If we use the simpler form with e = g = in Eq. (|T4*j) to fit the lowest 
three non-trivial moments, the parameters b and c that give the best fit to the data differ 
from those of the underlying distribution by approximately 10% and 30%, respectively. This 
provides a guide to the size of systematic errors associated with the choice of the parametric 
form. 

B. Valence distribution from lattice moments 

Having investigated the accuracy with which one can reliably extract the x dependence of 
the valence quark distribution from the lowest few moments, we now turn to the extrapolated 
lattice data discussed in Section lTTTl and use these to fit the parameters A, b and c in Eq. (fT^j) . 
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There are several possible approaches to reconstructing the x distribution from the avail- 
able data, which we discuss in the following. 

(i) Ideally, the n-odd and n-even moments should be fitted independently as they corre- 
spond to different distributions (Eq. ©), and the valence distribution extracted from the 
n-even (C-odd) lattice moments. In this approach (which we refer to as method (i)), both 
the statistical and systematic uncertainties (associated with the fact that current lattice 
data do not constrain /i) of the various extrapolations would be improved by future lattice 
data. However, the two available (n = 0, 2) moments are not sufficient to constrain all of the 
parameters in the standard form of Eq. (j!2j) . With the existing data, taking u = 0.7 GeV 
in the extrapolation, the n = moment fixes A, and we find a minimum x 2 along the line 
b ~ —0.9 + 0.2c (for < c < 4). For the case c = 1, one has b ~ —0.7, while for c = 2, 
b ~ —0.5. Both of these curves are in qualitative agreement with the Drell-Yan data. If it 
proved feasible to extract the n = 4 and 6 lattice moments in the future, this method would 
be ideal. 

(ii) An alternative approach is to assume, as discussed in Sec. that quark loops are 
suppressed at large m n and that the n-odd valence moments are approximately equal to the 
calculated C-even moments in that mass region. This provides us with 4 valence moments 
to which we fit 3 parameters. We choose fi = 0.7 GeV for the central extrapolation, taking 
H = 0.4 and 1.0 GeV and the lattice data ± its quoted errors, respectively, as a conservative 
measure of the overall error. This yields the parameters shown in Table ITT1 (method (ii)). We 
determine the uncertainty in the parameters (arising from the statistical errors in the lattice 
data and the systematic errors in the extrapolation (choice of fj,)) by choosing an ensemble of 
sets of randomly distributed moments within the extrapolated bounds and computing their 
standard deviation. As discussed above, there is additional systematic uncertainty arising 
from the reconstruction procedure that is not shown. 

The resulting x distribution, which is illustrated in Fig. |2j is in qualitative agreement with 

n 

the Drell-Yan pion structure function data |5J]. At intermediate values of x, the extracted 
distribution tends to lie a little above the experimental data, while at large x it appears 
to lie slightly below - with an x —>■ 1 behavior similar to that predicted by hadron helicity 
conservation in perturbative QCD. Of course, given the relatively large errors at present (the 
shaded region is the envelope of the ensemble of reconstructed distributions used in the error 
analysis), the distribution shows no significant disagreement with the experimental data. 
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FIG. 2: Valence distribution of the pion, reconstructed using method (ii) described in the text. 
The shaded region corresponds to the uncertainty in the distribution (see text). The dashed 
line corresponds to v n (x) = J 4x~ 1 / 2 (1 — x) 2 , which incorporates the hadron helicity conservation 
expectation for the large- x behavior. Experimental Drell-Yan data (diamonds) are from Ref. 

Since Eq. © gives the m n dependence of the moments, we can examine the dependence 
of the valence distribution on the pion mass. The result of reconstructing the PDF at 
several values of m 2 is shown in Fig. 01 We do not show results for values of m 2 larger 
than 1 GeV 2 because there are no lattice data to constrain the reconstruction. However, 
we have checked that if the heavy quark limit is built into our extrapolation function, the 
distribution approaches a S-f unction at x = 1/2. It is interesting that, even without such 
a constraint, the PDF seems to show that the momentum of the pion is shared primarily 
between the two valence quarks for above 0.7 GeV. 

The obvious problem with this method is that the assumption that the C-odd and C-even 
moments are approximately equal must break down as the lattice data are extended to lower 
masses — presumably where one begins to see curvature in the n = 2 moment. 

(iii) A third possibility is to extrapolate the n-odd moments linearly, according to Eq. (JJJ), 
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FIG. 3: Mass dependence of the valence quark distribution of the pion (with in GeV). 

and subtract twice the moments of the phenomenological sea at the physical mass to give 
the valence moments. The disadvantage of this method is that it relies on phenomenolog- 
ical information in addition to that obtained from the lattice. Moreover, the sea quark 
distribution in the pion is only very weakly constrained by experimental data, so that the 
errors on the valence, n-odd moments will be large compared to those on the n-odd total 
moments extracted from the lattice. A further disadvantage of this method, from the purely 
theoretical point of view, is that in this approach one obviously cannot study the variation 
of the pion PDFs as a function of quark mass, as was done for method (ii) above and for 
the nucleon in Ref. On the other hand, the errors on the extrapolations can be im- 



proved systematically as the lattice data at sma! 



ler quark masses become available and new 



experiments better constrain the pion sea P. Is3|. 

Taking the n-odd sea moments from an average of the SMRS and GRS distributions 
and extrapolating the n = 2 moment with [i = 0.7 GeV gives the parameters shown in 
Table ITT1 as method (iii). Errors are as described for method (ii) (where relevant). These 
parameters yield a distribution which is in good agreement with the Drell-Yan data, as seen 
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FIG. 4: As in Fig. |2J but using method (iii) to reconstruct the valence pion distribution. 

in Fig. |3J In particular, the extracted curves are consistent with (though slightly harder 
than) the x — > 1 behavior found in the experimental analyses, which, however, disagree with 
the hadron helicity conservation predictions. 

While there is some difference between the detailed x dependence for the valence quark 
distribution obtained using the methods (ii) and (iii), these should disappear once more 
accurate lattice data on the moments become available. Nevertheless, the fact that both 
methods are in reasonable agreement with the Drell-Yan data, and with each other, is very 
encouraging. It would be particularly valuable to have accurate higher moments (n = 4-6) 
in order to constrain the detailed shape of the distribution. 

C. Sea distribution 

The sea quark distribution in the pion, defined as 



(15) 
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Fit A b c 

SMRS 1.08 -0.36 1.08 

GRS 0.98 -0.47 1.02 

method (ii) 4.4 0.1(5) 2.5(1.5) 

method (iii) 0.6 -0.6(3) 0.8(9) 



TABLE II: Parameters of the valence distributions from various methods of analysis. The GRS 
valence distribution ()13|) in addition uses the parameters e = —0.81 and g = 0.64. The quoted errors 
combine the statistical and systematic (from fj,) errors. (Note that the lower limit on the parameter 
c for method (iii) is constrained to be positive.) Errors are not given for A as it is constrained by 
normalization. There is some additional systematic uncertainty from the reconstruction procedure 
which is not shown but is discussed in the text. 

is relatively poorly known experimentally. There are no data from the Drell-Yan reaction 
for x < 0.2 P, 0], and the size of the sea is constrained only by imposing the momentum 
sum rule. A simple parameterization of the pion sea (as used by SMRS |ll() is 

x s^x) = A s (l - x) v . (16) 

The analysis of GRS determines the pion sea with reference to the proton sea. A similar 
constraint can be derived at x ~ 10 -2 from semi-inclusive DIS measurements at HERA, 



with the result F? » F.f/3 |53(. This finding tends to favor the SMRS sea over that of 
fl 

GRS |6j , however, this information corresponds to such small values of x that it is of little 
assistance for the present analysis. 

In order to obtain information on the pion sea from the lattice data, one would ideally 
extract the valence distribution according to method (i) above, and use this to calculate the 
n-odd valence moments. These would then be subtracted from the total moments, obtained 
from the corresponding extrapolation of the lattice data to obtain the n-odd moments of 
the pion sea at the physical pion mass. The x dependence of the sea distribution could then 
be reconstructed using the form (|16|). given enough moments. 

In the absence of sufficiently many moments for this to be a practical solution, an alterna- 
tive is to use the phenomenological valence (n-odd) moments instead of the lattice moments. 
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Since these moments are relatively well known, this procedure should be reasonably reliable. 
The n-odd sea moments which we extract using the linearly extrapolated total moments 
from the lattice minus the phenomenological valence moments are given in Table HI Using 
these data to fit the Mellin transform of Eq. (JTEJ), we find the parameters A s = 0.27 and 
f] = 5.8. Unfortunately, the statistical uncertainty in these lattice sea moments is large and 
the constraints on these parameters are very weak. In particular, the third moment of the 
sea is consistent with zero: for (a; 3 ) sea — > the reconstruction gives rj — ► oo. Nevertheless, 
this is in principle improvable and if the size of the errors on the n = 3 lattice moment were 
comparable to that on the n = 1 moment, a more robust reconstruction could be performed. 

One could also modify this method by including the n = 2 sea moment constructed 
from the difference of the linear and chiral extrapolations of the lattice data. However, 
this introduces additional uncertainty (from //) in the analysis and does not reduce the 
uncertainty in the reconstructions. 

V. SUMMARY AND PROSPECTS 

We have studied the problem of the chiral extrapolation of the moments of the PDFs of 
the pion, from the large quark masses where current lattice QCD calculations are performed 
to the physical values. As in earlier studies of the PDFs of the nucleon, the non-linearity 
of the model independent non-analytic variation of the moments of the valence PDF is 
extremely important, producing a significant change in the moments at the physical quark 
mass, compared with a naive linear extrapolation. In comparison, the moments of the singlet 
distribution, q + q, show no non-analytic behavior and are therefore expected to extrapolate 
smoothly to the chiral limit. 

Having studied the extrapolation of the moments of both the singlet and non-singlet 
PDFs, we examine various procedures for reconstructing the valence and sea distributions 
of the pion from the lattice moments. To make optimal use of the available lattice data 
for the n = 1-3 moments, we make the reasonable assumption that at large quark masses 
(m % ^ 500-600 MeV) the effects of quark loops are suppressed, so that the effects of the 
quenched approximation and disconnected insertions will not affect the extrapolation. This 
allows the parameters of the valence, and to some extent the sea, quark distributions to 
be determined, and the extracted distribution compared with the available Drell-Yan data. 
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Over the range of intermediate x, from 0.2 to 0.8, the reconstructed valence distribution is 
in fair agreement with the Drell-Yan data, within the rather large errors arising from the 
extrapolation procedure. At this stage, however, it is not possible to distinguish between 
the large- x behavior predicted from hadron helicity conservation in perturbative QCD and 
that found in the Drell-Yan data. Nevertheless, new lattice simulations, on the much faster 
computers expected to be devoted to lattice QCD in the next few years, should offer the 
chance, when analyzed using the techniques set out here, to determine leading twist PDFs 
with an accuracy that exceeds that of current experiments. 

The results of this analysis can be used to guide future studies of PDFs in lattice gauge 
theory. Specifically, 

• The clearest observation is that it would be extremely valuable to have quenched 
calculations for the n = 2 and n = 3 moments of an accuracy comparable to that for 
n = 1. This is especially important for pion masses below 0.4 GeV 2 . This alone would 
make a substantial improvement in the errors on the parton distribution functions. 
To better constrain the functional form of the ^-dependence, calculations of several 
higher moments (e.g. n = 4,5) would also be desirable. 

• In the case of the nucleon there is no observable difference between the moments 
calculated in quenched and full QCD in the mass range covered. It is vital to check 
that this is also true for the pion. 

• In order to better constrain the extrapolation and to assure us that we are on the right 
track it is important to push the lattice simulations to lower pion masses. Ideally this 
should occur in full (unquenched) QCD; however, quenched, and especially partially- 
quenched (which is a computationally efficient way to get to smaller valence quark 
mass without actually ignoring quark loops) simulations would also provide valuable 
information to guide the chiral extrapolation. 

• Of course, if one wants to explore the sea quark distribution, for which there is little 
hard information at present, it will also be necessary to include disconnected quark 
loops, even though it is difficult to pick out a signal for such terms J3|. 

• Finally, as with all lattice simulations, we need to confirm that the continuum (a — > 0) 
and infinite volume (L — > oo) corrections are fully under control. 
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With such a program we could expect significant advances over the next few years in our 
understanding of the quark structure of the pion in QCD. 
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